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1 Introduction 



D°D° oscillations arise from beating of two mass eigenstates, Di, D 2 that propagate 
differently in time, and are related to the instaneous flavour states by 



|£>i) = p\D°} + q\D ) 
\D 2 ) = p\D°}-q\D°) 



| bl 2 + \l\ 2 = 1 ; arg{g/p} = M . (1) 



Their amplitude and frequency depend upon the normalized differences 

x = (m 1 -m 2 )/V ; y = (T 1 - T 2 )/(2T) ; V = (I\ + T 2 )/2, (2) 

in the D\ and D 2 masses and decay rates, and Ti^, respectively. 

The leading term in the standard model, SM, contributing to AC = 2 transitions 
(C is the charm quantum number) is a box diagram with W and d-type quarks. This 
predicts very little mixing (x < 1CT 5 ) pQ, because the SM was designed [2] to produce 
large cancellation, ~ (m 2 s — mfy/ml, between s and d quarks, and also because the 
Cabibbo-Kobayashi-Moskawa (CKM) coupling V u b is small. New physics (NP) could 
lead to larger values for x or y |3j SI EJ , and could also produce CP violation {CPV) in 
mixing. However, long range SM effects from real, AC = 1 intermediate hadron states 
have also been suggested [B] to contribute significantly to both x and y. Theoretical 
uncertainty in ways to estimate these sums exists [7], and values for x and y over a 
wide range (10 _7 -10 -2 ) have been quoted [3l|8]. The present experimental results fit 
well into the higher range, but the possibility that NP may be involved cannot be ruled 
out. Theoretical consensus is that NP is neither required nor ruled out by present 
measurements. It is mostly agreed that observable CPV at present experimental 
sensitivities would signify unambiguous evidence for NP were it observed jl]. 

1.1 D° decays 

Measurements of D°D° mixing have been based on observed time-dependences for 
decays of D° to final states / accessible to either D° or D° (e.g. K + K~, K + , K + 
7T~ 7T°, etcQ. In such decays, mixing and direct decay interfere. In the abscence of 
CPV, the number N(t) of D° remaining at time t is given, to second order in x and 
V by 



N(t) = N(0)e- Tt x 



2 2 

" - l + lA^^cos^-xsin^^ + ^^IA;! 2 ^) 2 



(3) 



The first term in square parentheses corresponds to direct decay D° — > f (amplitude 
A/). The term quadratic in t represents mixing (D° — > D°) followed by decay D° — > f 



* Unless explicitly stated otherwise, charge conjugate states are implied. 
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(amplitude Af). The middle term, linear in t, is due to the interference between these 
processes. The parameter Xf is defined as 

\ f = (qA f ) I (pA f ) ; arg{A/} = 9 f = <f> M + (j) f + S f . (4) 

These amplitudes Af and Af have relative weak (strong) phases 0/ (5f). 

If CP is conserved in mixing, p = q = l/v2 and <p M = 0. If it is conserved in 
decay, then (fif = 0. In the special case that / is a CP-eigenstate, then (fif = Sf = 
and Xf = ±e l< f >M . In most measurements so far, CPV is ignored and then Of = 5f. 

The strong phase difference 5f is, in general, unknown and means that, for many 
analyses, only the rotated parameters 

x = x cos 5/ + ysin5f ; y = ycos5f — xsin5f (5) 

can be determined. 



1.2 Flavor tagging and decay lengths 

Identifying the flavor (D° or D°) of a neutral D at birth is sometimes required. When 
this is so, each D° candidate is required to come from a D* + — > -D°7r+ decay, where 
the sign of the low momentum pion, ir s , tags the D flavor. D*'s are identified by the 
required peak in the distribution of AM, the difference between the invariant mass M 
of the D° daughters and that of the -D°vr+ system. This peak also serves to improve 
on signal selection and background rejection. 

Decay length distributions for D 0, s from B decays differ from those produced 
directly. In the Belle and BABAR analyses, D 0, s from B are removed by selecting 
center of mass momentum of the D° above the kinematic limit (~2.5 GeV/c). In 
the CDF data, the two sets of decays are distinguished by the distributions of their 
impact parameters. 



2 Evidence for mixing in decays to CP eigenstates 

The Belle collaboration first reported evidence for mixing with 3.2a significance in 
decays of D° to CP eigenstates h~h + {h = n or K) [9]. For such decays, if CP 
is conserved, mean decay times computed from Eq. ^ r are related to decays to 
non-CP states such as D° — » K~n + f\by 

t(D° K-TT+) 



^Here, as in all cases unless stated otherwise, the first uncertainty is statistical and the second is 
systematic. 
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Using separate, flavor-tagged samples of D° and D°, Belle also measured the CP 
asymmetry: 

= x( C ° -> ^-) - -> ^-) = 

consistent with zero. 

This evidence was confirmed by the BABL4R collaboration in two independent mea- 
surements using flavor-tagged D° decays to K~^K~ and 7r + 7r~ [TU] and to a large, 
disjoint, untagged K + K~ sample [H] with a combined significance of 4.1er. The CP 
asymmetry A T for the BABAR tagged decays was 0.260±0.360±0.080%, also consistent 
with zero. 

These results, with other less significant ones have been combined [12] to give a 
mean value yep = (1.107 ± 0.217) x 10~ 2 , 5. Oct from the no mixing value zero. 



3 Evidence for mixing in "wrong sign" decays 

"Wrong sign" (WS) decays D° — > K + 7i^(n°) are sensitive to mixing. Here, direct 
decay is doubly Cabibbo suppressed (DCS) and the three terms in Eq. ^ are all 
comparable. For Cabibbo-favoured (CF) "right-sign" (RS) decays D° — > K~ n + (7r°) , 
however, the first term dominates and decays have an exponential distribution. Ex- 
perimentally, decay time distributions for both rare WS and copious RS events are 
fit simultaneously to determine the time resolution parameters (same for both) and 
mixing parameters (only observable in the WS sample). For D° — > K + n~ decays, Sf 
is unknown and only rotated parameters x' 2 and y' can be measured. 

The BABAR collaboration was first to report evidence for D° oscillations [13] in 
a large sample of WS decays D° — > K + tt~ . In a challenging and careful analysis 
of the decay time distribution in which the time resolution was only slightly less 
than the oscillation period, they reported a 3.9a deviation. In an earlier analysis of 
this channel, Belle [H] had seen only a 2. Oct effect [^J Subsequently, the evidence for 
mixing was confirmed by CDF [15] with a 3.8<r significance. Results from these three 
experiments are summarized in Table [T] 



3.1 Mixing in D° — > K + ir~ir° decays 

The BABAR collaboration [16] also reported evidence for mixing from a time-dependent 
Dalitz plot (DP) analysis of WS decays to the three-body system having an additional 
7T° meson. This system is similar to the two-body one, except that the final state / 

■•■This analysis was even more sensitive than that of BAR4R, but the central values for x' 2 and y' 
were closer to zero. 
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Table 1: Rotated mixing parameters (x' 2 ,y') in units of 10 -3 from fits to WS D° — > 
K + 7i~ decays described in the text. The significance for mixing is obtained from 2- 
dimensional likelihood contours between the highly correlated x' 2 and y' values, that 
are computed to include systematic effects. The significance is taken from the contour 
on which the no mixing point (x' 2 = y' = 0) lies. The table includes CP asymmetries 
for the DCS/CF ratio R D = |A|~ 2 for D° and R D = |A|~ 2 for D°, and mixing rates 
R M = x' 2 + y' 2 for D° and R M = x' 2 + y' 2 for D°. 

Parameter BABAR CDF Belle 

x' 2 -0.22 ±0.37 -0.12 ±0.35 0.18+^ 

y' 9.7 ±5.4 8.5 ±7.6 0.6^ 

Mixing significance 3.9a 3.8a 2.0a 

a D = (R D -R D )/(R D + R D ) -21 ±54 23 ± 47 

a M = (R M - R M )/{R M + R M ) 670 ± 1200 



is now a point in the K + ir~ir° phase space, specified by its DP coordinates So and s + , 
the squares of invariant masses of K + ir~ and K + ir° systems, respectively. 

The expected time-dependence is also given by Eq. (|3]), with the parameter A/ 
and its phase at each / given by the ratio of decay amplitudes A(sq, s+) for CF decay 
(D° K+n-n ) and A(s , s+) for direct DCS decay (D° -»■ K + 7T~7r°) at that point. 

In the BABAR analysis, approximately 660K RS (R-n+n ) and 3K WS (K + n-n°) 
flavor-tagged decays are extracted from a 384 ftT 1 sample of e + e~ interactions. Pa- 
rameters for separate models Af and Af for CF and DCS amplitudes, respectively, in 
each case a linear combination of Breit-Wigner (BW) amplitudes describing their dif- 
ferent Kir and irir resonance structures, are determined from simultaneous fits to the 
DP's for these samples^ Assuming CP is conserved, a model for the time-dependent 
variations in 5f from point to point in the DP are thus obtained, but an unknown 
overall phase Sk-k-k allows only rotated coordinates x" and y", similar to x' and y' in 
Eq. (pjl), to be measured. 

Values obtained for x" and y" are given in Table [2j A confidence level test, similar 
to that used in the WS K + tc~ case indicates evidence for mixing at the 3. la level. 
The major systematic uncertainties are associated with the assumptions in the BW 
model and in the description of the large background under the WS events. Estimates 
of these are included in the computation of the confidence level. 

The fit procedure is repeated separately for D° and D° samples to obtain values 
for x" and y" also listed in Table |2j indicating no evidence for CPV . 

§ Since the RS sample is dominated by CF decays, a time- integrated fit is made to determine Af. 
A time-dependent fit to the WS sample determines parameters for Af and x" and y" values. 
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Table 2: Rotated mixing parameters x" and y" from fits to BABAR data described in 
the text. The first error is statistical and the second that attributed to systematic 
effects. 

Sample x" (%) y" (%) 

D° and D° 2.61^ ± 0.39 -0.06lg;g ± 0.34 
P>° only 2.53^ ± 0.39 -0.05lg;g ± 0.50 
/J° only 3.55±g;^ ± 0.65 -0.54+°;^° ± 0.41 



4 Decays £>° K^ttt^ and £>° ifOif+if- 

Decays to self-conjugate systems (sum of CP-even and CP-odd states) allow direct 
measurement of x and y since the strong phase difference between D° and D° decays is 
zero. Using these channels, it is also possible to obtain values for the CFV parameters 
\q/p\ and <p M . 

This was first exploited for the D° — > K®tt + it~ mode by CLEO [T7j using a 9 fb _1 
data sample and obtaining only an upper limit on mixing. The Belle collaboration 
repeated the analysis [IB] with a 540 fb _1 sample to obtain central values for x, 
y, \q/p\ and arg{g/p}. More recently, BABAR has used both D° — > K®7m~ and 
D° K®K + K~ modes for a 486.5 fb _1 sample to obtain the most precise results. 
Parameters from all three experiments are summarized in Table |3j 



Table 3: Mixing parameters from fits to D° — > K Q s h + h~ (h = tt or K) decays. The first 
uncertainty is statistical, the second is from systematic effects. A third uncertainty 
comes from ambiguities in the choice of model used to describe the decay amplitudes. 



Experiment 




Mixing 


Parameters 


CLEO 2.5 (9 fb^ 1 ) 


X- 


=( 1.9+^ ±0.4 ±0.4)% 


y=(-1.4±2.4±0.8±0.4)% 


BELLE (540 fb" 1 ) 
(Allowing CFV) 


X- 

\q/p\-- 


= (0.81 ±0.30^-1^)% 
=0-86185 Tot ± 0-08 


y=(0.37±0.25i[J- 


BABAR (486.5 ftT 1 ) 


X- 


= (0.16 ±0.23 ±0.12 ±0.08 


)% y=(0.57± 0.20 ±0.13 ±0.07)% 


(D° only) 


x + - 


= (0.00 ±0.33)% 


y+=(0.55±0.27)% 


(D° only) 


x~- 


= (0.33 ±0.33)% 


y-=(0.59 ± 0.28)% 



These results agree well and, since they do not depend upon any unknown strong 
phases, significantly affect the averages. Central values for x and y, however, are such 
that the significance for mixing is small (2.2<r for Belle and 1.9<r for BABAR). 

The decay amplitude models describing the DP distributions differ between the 
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three K®ii + ii~ analyses. We note from Table |3| that uncertainties from the decay 
amplitude model observed in BABAR and Belle lead to irreducible uncertainties of 
order 10~ 3 in both x and y. These, in turn, limit CPV measurements to ~25% in 
\q/p\ - 1. 

5 Future Outlook 

The task of combining some 30 "mixing observables" , some of which are presented 
above, into values for the important physics parameters underlying them has been 
undertaken by the heavy flavor averaging group (HFAG [12] ). These are summarized 
in Table [4| \ 2 contours for this fit indicate that these central values are far (at least 
lOcx) from the no mixing point (x — y — 0) and are within la from the no CPV point 
(|g/p| = l,0 M = O). 

Table 4: HFAG summary of mixing parameters from fits to 30 observables. 



X 


(6.3t^) x 10- 3 


y 


(7.5 ± 1.2) x 10~ 3 


\q/p\ 


Q1+ - 18 

u - yi -0.16 


fit 


-10 2+ 9 - 4 


R D 


(3.309 ±0.081) x 10~ 3 


A D 


(-19.2 ±24) x 10~ 3 




22.01& 


K 

K.7T1Y 


IQ Q+21-8 
iy.o_ 22 .9 



Upgrades in both B factories are planned that will increase event yields by a factor 
~ 100. New results from D°D° pair production at charm threshold from BES III are 
also anticipated to provide improved measurements of strong phase differences Sf for 
K + 7i~ , K + 7t~7t° and also for individual points in various Dalitz plots that could be 
used for reducing uncertainties in decay amplitude models. LHCb is running now, 
too, and will surely add to our knowledge and precision of mixing parameters very 
soon. 

As we look ahead, we can consider how we might search for CPV, a possible 
indicator for NP. One way is to use the self-conjugate channels D° — > K®h + h~ to 
directly measure \q/p\ and 0a/- From Table [3] it is seen that these are limited by 
decay amplitude model uncertainty to ~ 8% and 3°, respectively. Another way that 
has been considered in detail by the Superl? collaboration [19] is to measure CP 
asymmetries in the parameters x, y, x', y', x" or y" that each provide measurement 
of \q/p\ 2 — 1. To reach a precision in \q/p\ of 1% requires precision in x or y of ~ 10~ 4 
that can be acheived only if decay model uncertainties can be reduced by a factor 
10, as illustrated in Fig. [TJ Asymmetries in these parameters in different channels at 
the 5% level are also possible, and can be used to test for possible presence of direct 
CPV. 
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Figure 1: Mixing observables from BABAR projected into the (x,y) plane. Shaded 
areas indicate the 68.3% confidence region for each. Contours enclosing 1 — 5a two- 
dimensional confidence regions from x 2 fits to these observables are drawn as solid 
blue lines, (a) Includes results anticipated from scaling BABAR results by a factor 
100 in statistical significance anticipated for Super5. Uncertainties in x(±7.5 x 10~ 4 ) 
and j/(±1.9 x 10~ 4 ), are dominated by uncertainty in the amplitude model. In (b), 
projections of hypothetical measurements of Sk+w- and Sk+^-^o are included and 
are expected to reduce the amplitude model uncertainty by a factor 10 leading to 
uncertainties in x(±2.0 x 10" 4 ) and y(±1.2 x 10~ 4 ). 



6 Summary 

In conclusion, there is strong evidence for mixing from a variety of measurements, 
but no one observation at the 5a level has yet been made. No evidence has yet been 
seen for CFV . Prospects are that Superl? or Belle II will improve precision in x and 
y to the 10~ 4 level providing sensitivity to CFV (\q/p\) at about 2 — 5% from channel 
to channel. However, this will only be possible if either new measurements of strong 
phases, or new amplitude analysis techniques become available. 
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